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Abstract 

We present an updated version of the SimProp Monte Carlo code to study the propagation of ultra high energy 
cosmic rays in astrophysical backgrounds computing the cosmogenic neutrino fluxes expected on earth. The study of 
secondary neutrinos provides a powerful tool to constrain the source models of these extremely energetic particles. We 
will show how the newly detected IceCube neutrino events at PeV energies together with the the latest experimental 
results of the Pierre Auger Observatory and Telescope Array experiment are almost at the level of excluding several 
hypothesis on the astrophysical sources of ultra high energy cosmic rays. Results presented here can be also used to 
evaluate the discovery capabilities of future high energy cosmic rays and neutrino detectors. 
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1. Introduction 

The study of Ultra High Energy Cosmic Rays 
(UHECR) started already in ’60 with the first pioneer¬ 
ing observations of the Volcano Ranch experiment, that, 
in 1962, observed the first cosmic ray event with energy 
larger than 10 2l) eV (TJ. Nowadays the most evolute ex¬ 
periments observing UHECR are the Pierre Auger Ob¬ 
servatory in Argentina 12, far the largest experimen¬ 
tal setup devoted to the study of these particles, and 
the Telescope Array (TA) experiment |3), placed in the 
United States, with roughly 1/10 of the Auger statistics 
at the highest energies. 

The experimental study of UHECR clarified few im¬ 
portant characteristics of these particles: (i) UHECR 
are charged particles with a limit on photon and neu¬ 
trino fluxes around 10 19 eV at the level of few percent 
and well below respectively aim (ii) the spectra 
observed on earth show a slight flattening at energies 
around 5 x 10 18 eV (called the ankle) with (iii) a steep 


suppression at the highest energies - 10 2 ° eV 1213. 

The propagation of UHECR from the source to the 
observer is conditioned by the expansion of the universe 
and the interaction with astrophysical backgrounds, 
namely the Cosmic Microwave Background (CMB) and 
the Extragalactic Background Light (EBL). While the 
propagation of nucleonf] is conditioned only by the 
CMB field the propagation of nuclei is also affected 
by the EBL. Apart from the expansion of the universe 
that, adiabatically, reduces the energy of any propagat¬ 
ing particle, the processes that involve protons are: (i) 
pah production and (ii) photo-pion production; those in¬ 
volving nuclei are: (i) pair production and (ii) photo¬ 
disintegration |8][9]]. 

One of the most important observables in the physics 
of UHECR is certainly the chemical composition of 


'Hereafter discussing freely propagating UHE nucleons we will 
always refer only to protons because the decay time of neutrons is 
much shorter than all other time scales involved GD0- 
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these particles. The experimental observations of the 
composition are not conclusive, with different results 
claimed by TA and Auger. The analysis performed by 
TA is compatible with a proton dominated composition 
at all energies, starting from the lowest around 10 18 eV 
up to the highest at 10 20 eV |3) . On the other hand the 
observations of Auger show a more rich phenomenol¬ 
ogy with the lowest energies dominated by protons and, 
starting from energies around 3 x 10 18 eV, a composi¬ 
tion more and more dominated by heavier nuclei with 
a strongly reduced number of protons at energies above 
2 x 10 19 eV J2|. 

Restricting the analysis to a pure proton composi¬ 
tion, as appropriate to the interpretation of TA data, 
the UHECR observations can be elegantly explained 
by the dip model IflOl 11 11. In the framework of this 
model the flux behaviour in the region of the ankle 
is due to the effect of the proton pair-production pro¬ 
cess on the CMB radiation field iflOl . While the strong 
flux suppression at the highest energies is the effect 
of the photo-pion production process, the so-called 
Greisen-Zatsepin-Kuzmin (GZK) cut-off HI 21 11 311. Tak¬ 
ing into account the sole Auger data, to reach a rea¬ 
sonable agreement with observations, one should con¬ 
sider mixed compositions. Assuming that protons give 
their principal contribution only at the lowest energies 
< 3xio 19 ev mm, below the photo-pion production 
threshold ~ 6 x 10 19 eV. In this case the flux suppres¬ 
sion observed at the highest energies is due to the photo¬ 
disintegration process suffered by heavy nuclei mm. 

The production of secondary particles due to the 
interactions of UHECR with background photons is 
strongly dependent on the chemical composition. As 
was first realised in ifTTl . the proton content at the high¬ 
est energies is a crucial quantity that regulates the fluxes 
of secondary (cosmogenic) photons and neutrinos. In 
the present paper, implementing an updated version of 
the 5/mPro/^jMonte Carlo (MC) code jl 81 , we will dis¬ 
cuss the production of secondary neutrinos by the prop¬ 
agation of UHECR. The study presented here, discussed 
in detail in m, has a twofold interest: from one side 
using the latest observations of the IceCube Il20l and 
Auger 0 we can already draw interesting conclusions 
on the sources of UHECR; on the other side this study 
should be intended as a benchmark computation to asses 
the discovery capabilities of the next generation experi¬ 
ments. 


2 SimProp is available upon request writing to the authors or at 
SimProp-dev@ aquila. infn. it. 


2. UHECR models and secondary neutrinos 

There are two processes by which neutrinos can be 
produced in the propagation of UHECRs: (i) the de¬ 
cay of charged pions produced by photo-pion produc¬ 
tion, 7^ —> p ± 4- and the subsequent muon de¬ 

cay yu* —» e ± + v^fvyi) + v e (v e ); (ii) the beta decay of 
neutrons and nuclei produced by photo-disintegration: 
n —> p + e~ + v e , (A,Z) —» (A.Z - 1) + e + + v e , or 
(A.Z) —> (A,Z + 1) + e~ + v e . These processes pro¬ 
duce neutrinos in different energy ranges: in the former 
the energy of each neutrino is around a few per cent 
that of the parent nucleon, whereas in the latter it is less 
than one part per thousand (in the case of neutron de¬ 
cay, larger for certain unstable nuclei). This means that 
in the interactions with CMB photons, with a Lorentz 
factor threshold around T > 10 10 , neutrinos are pro¬ 
duced with energies of the order of 10 18 eV and 10 16 eV 
respectively. Interactions with EBL photons contribute 
with a much lower probability respect to CMB photons, 
affecting a small fraction of the propagating protons and 
nuclei. Neutrinos produced through interactions with 
EBL, characterised by lower thresholds, have energies 
of the order of 10 14 eV in the case of photo-pion pro¬ 
duction and 10 16 eV in the case of neutron decay. 

Neutrinos produced by the interaction of UHECR, 
because of their extremely low interaction rate, arrive 
on Earth unmodified with the overall universe contribut¬ 
ing to their flux. This is an important point that makes 
neutrinos a viable probe not only of the chemical com¬ 
position of UHECR but also of the cosmological evolu¬ 
tion of sources that, as we will show below, can be also 
constrained by the neutrino flux observed on Earth. 

In the following we will consider the two cases of the 
dip model and mixed composition, discussing the ex¬ 
pected neutrino flux with different assumptions on the 
cosmological evolution of sources. We will consider 
the case of sources (see mi and reference therein) (i) 
with no cosmological evolution, (ii) with the same cos¬ 
mological evolution of Active Galactic Nuclei (AGN), 
supposed to play a role in particles acceleration till the 
highest energies m, and (iii) with the cosmological 
evolution of the Star Formation Rate (SFR). All com¬ 
putations presented here, discussed in detail in m . are 
performed under the assumption of a homogenous dis¬ 
tribution of sources. This assumption does not affect the 
expected neutrino spectra because in the case of neu¬ 
trinos the overall universe, till the maximum red-shift, 
contributes to the fluxes. Possible flux variations due 
to a local inhomogeneity in sources distribution gives a 
negligible contribution to the total flux. We also fix a 
maximum red-shift of the sources z max = 10, which is 
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Figure 1: Fluxes of UHECR expected at Earth in the case of the dip 
model (pure protons). Fluxes are computed in the three cases: AGN 
cosmological evolution y g = 2.4 dot-dashed line, SFR evolution with 
y g = 2.5 dashed line and no cosmological evolution with y g = 2.6 
solid line. Experimental data are those of TA (3| and Auger (2)- 

the typical red-shift of the first stars (pop III) EH . in 
any case if z max > 3 the expected fluxes of primary and 
secondary particles is almost independent of z m ax Ifl9l . 

Once produced at cosmological distances neutrinos 
travel toward the observer almost freely, the opacity of 
the universe to neutrinos starts to be relevant only at the 
highest red-shifts namely at z > 10 Il22ll23l . Therefore, 
given the assumptions discussed above, in our computa¬ 
tions we have neglected any effect due to neutrino prop¬ 
agation apart from the adiabatic energy losses due to the 
expansion of the universe mm. 

As discussed in detail in m, the case of the dip 
model gives the highest flux of secondary neutrinos as it 
follows from the high efficiency of the photo-pion pro¬ 
duction process that is the major responsible for the pro¬ 
duction of secondary particles. In order to reproduce the 
TA observations the injection power law index y g needs 
to be changed according to the cosmological evolution: 
we used y g — 2.6, 2.5, and 2.4 respectively in the three 
cases of no cosmological evolution, SFR and AGN. The 
result in terms of UHECR flux is plotted in figure |T| 
together with the TA data. In figure [2] we plot the corre¬ 
sponding neutrino flux in the three cases of cosmolog¬ 
ical evolution: no evolution (red band), SFR evolution 
(green band) and AGN evolution (blu band). The flux of 
neutrinos depends on the assumptions made about the 
EBL cosmological evolution, therefore we have plot¬ 
ted the corresponding fluxes as coloured band to high¬ 
light the uncertainties related to the EBL evolution, the 
EBL evolution models used are those of Stecker m 
and Kneiske l25l . 

Results presented in figure [2] are quite interesting. 



Figure 2: Fluxes of neutrinos expected at Earth in the dip model 
with the three different assumptions on the cosmological evolution 
of sources, from top to bottom: AGN evolution blu band, SFR evolu¬ 
tion green band, no evolution red band. Fluxes are plotted as coloured 
bands to highlight the uncertainties connected with the EBL cosmo¬ 
logical evolution assumed (see text). Experimental data on neutrino 
fluxes are by IceCube ED and Auger (6) as labeled. 


showing that models with cosmological evolution of 
sources stronger than the case of AGN predict a total 
neutrino flux (blu band in figure[2]) almost at the level of 
the experimental limits fixed by Auger and the observa¬ 
tions of IceCube. Moreover, cosmogenic neutrinos can 
explain the IceCube events only at PeV energies not be¬ 
low. 

Let us now consider the case of a mixed composi¬ 
tion, as appropriate to the interpretation of the Auger 
data. As shown by 1(141 scenarios with mixed compo¬ 
sition have peculiar characteristics that, in order to de¬ 
scribe Auger observations, imply different source fam¬ 
ilies. At the highest energies (E > 5 x 10 18 eV) Auger 
observes predominantly heavy elements, with A > 12, 
while at the lowest energies, as a matter of fact, all 
UHECR experiments observe a light composition dom¬ 
inated by protons. To reproduce such observations two 
classes of sources are needed: one that injects only light 
elements with low maximum energy and steep injection 
y g > 2.3 and one that injects also heavy elements with 
higher maximum energy and an extremely flat injection 
y g = \m. In the case of nuclei cosmological evo¬ 
lution severely constrains sources injecting heavy ele¬ 
ments. In this case the effect of photo-disintegration 
of nuclei produces too many secondary protons already 
assuming a cosmological evolution as for the SFR and 
AGN cases. Therefore, in this case, we have considered 
only cosmological evolution for sources providing the 
light component (p and He) leaving without evolution 
those sources providing heavy components. Results are 
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Figure 3: Fluxes of UHECR expected at Earth in the case of the mixed 
composition model of (H Experimental data are those of Auger m 
and TA (3- Cosmological evolution of sources is taken into account 
only in the case of light component with: AGN cosmological evolu¬ 
tion y g = 2.5 dot-dashed lines, SFR evolution with y g = 2.6 dashed 
lines and no evolution with y g = 2.7 solid lines. 

presented in figure [3] where we plot UHECR fluxes to¬ 
gether with Auger data, and in figure [4] where we plot 
the corresponding neutrino fluxes. Light elements are 
injected with a power law index y g = 2.7,2.6,2.5 in the 
three cases of no evolution, SFR and AGN evolutions 
respectively. As discussed in l T4lfl9l . this combination 
of sources and injections reproduces fairly well Auger 
observations on both flux (see figure [3]i and chemical 
composition cued, giving a neutrino flux which is 
substantially below the Auger observations and almost 
at the level of the IceCube observations at PeV energies 
in the case of AGN cosmological evolution and well be¬ 
low in the other cases of cosmological evolution. 

3. Conclusions 

We conclude stressing the importance of the obser¬ 
vation of secondary neutrinos produced by the propaga¬ 
tion of UHECR, as discussed in 02) the observation of 
secondary neutrinos can be extremely important in or¬ 
der to characterise the possible sources of UHECR. As 
follows from figure [2] in the case of the dip model, as¬ 
suming the AGN cosmological evolution, the neutrino 
flux is almost at the level of being detected by both 
Auger and IceCube. In the case of source models with 
acceleration of heavy elements, the production of sec¬ 
ondary neutrinos is quite suppressed respect to the case 
of pure protons, well below the detection thresholds of 
Auger and IceCube, this makes the non-observation of 
secondary neutrinos an indirect confirmation of the pic¬ 
ture emerging from Auger data. 


Figure 4: Fluxes of neutrinos expected at Earth in the case of the 
mixed composition model of m. The coloured bands refers to the 
three cases of cosmological evolution as discussed in the text, with the 
same color codes of figure [2] 
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